Entanglement as a quantum order parameter 
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We show that the quantum order parameters (QOP) associated with the transitions between a nor- 
mal conductor and a superconductor in the BCS and r;-pairing models and between a Mott-insulator 
and a superfluid in the Bose-Hubbard model are directly related to the amount of entanglement ex- 
istent in the ground state of each system. This connection gives a physical meaningful interpretation 
to these QOP, which shows the intrinsically quantum nature of the phase transitions considered. 
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It has recently become clear that entanglement plays 
an important role in the understanding of critical and 
thermodynamical properties of quantum systems 0, |3, 
Hi Hi IE IE Hi' The theory of entanglement, developed 
in the context of quantum information processing, has 
been applied in the study of properties of several im- 

Sirtant condensed matter systems, such as spin chains 
0, H S H IE ^'^d bosonic systems described by 
quadratic Hamiltonians 0]. Particularly, it was shown 
that near a quantum phase transition entanglement can 
be classified in the framework of scaling theory P, Q 
and that, using the so-called localizable entanglement, it 
is possible to define the concept of entanglement length 
which, in some cases, can detect phase transitions unde- 
tected by the traditional correlation length |^. In this 
letter we add another element to this picture, studying 
how entanglement is related to the order parameter (a 
quantity that assumes non-zero values in one phase, while 
is strictly zero in the other) in some important quantum 
models. 

Classical order parameters usually have a clear physi- 
cal interpretation, e.g., the density in the transition be- 
tween a liquid and a gas and the magnetization in the 
transition of a paramagnet into a ferromagnet. On the 
other hand, quantum order parameters (QOP) are not 
normally associated with any physical meaningful quan- 
tity. We show that the order parameters associated with 
the BCS model for normal superconductivity, with the 
77 pairing model for hight temperature superconductivity 
and with the Bose-Hubbard model for superfluidity are 
directly related to the amount of entanglement presented 
in those systems. This intriguing connection indicates a 
physical interpretation for (some) QOP: they quantify 
how quantum correlated are the parts of the system. 

The BCS Model: The BCS model is the most success- 
ful microscopic theory to describe superconductivity at 
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temperatures near to zero. Some properties of entangle- 
ment in the BCS model were studied in Refs. [Tol llll|. 
In Ref. [m in special, it was shown that entanglement 
is directly related to the order parameter in the reduced 
BCS model. We show that the this equivalence is in fact 
rather general for fermionic systems whose ground state 
is described by the BCS state. 

Following Ref. consider a fermionic system de- 

scribed by creation operators aj^, where the label a rep- 
resents the quantum numbers a = (fc,f7) of an electron 
and —a stands for the time-reversed state. The most gen- 
eral Hamiltonian for this type of system can be written 
as 

H^^Tapaiap + ^ ^ {a(i\V\-i5)a)^a!^pa^a5, (1) 

where Taj3 is the matrix element associated with the ki- 
netical energy and any external potential, T = P'^ /2M + 
W, and (a/3|y|7(5) is the matrix element from the anti- 
symmetric interaction ■ The crucial point in the BCS 
theory is the construction of an Ansatz for the funda- 
mental state of the system which takes into account the 
formation of Cooper pairs: 

\BCS) = l[{u^ + v^aialj\0), (2) 

where the real parameters Mq e Va satisfy 

ul+vl = l. (3) 

Is clearly seen that the BCS state consists of correlated 
electrons pairs, since the electrons described by a and 
—a are always associated. However, if either = 1 
or = 1 for every a, the ground state becomes sep- 
arable, i.e., the BCS state reduces to the Haartree-Fock 
(HF) approximation. The criterion for superconductivity 
is exactly the existence of a BCS state with lower energy 
than the HF state. Therefore, in a qualitative way, en- 
tanglement is necessary for superconductivity. 
It can be shown that the energy of the system is given 
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where n is the chemical potential, 



is the energy of each Cooper pair and 
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2[(e„-/i)2 + A^]i/2 



(4) 
(5) 

(6) 



The first term in the R.H.S. of equation (4) is just the HF 
energy. We thus see that the condition for the existence 
of superconductivity is that at least one of the is 
non-zero. The A^ are considered the order parameters 
and have quantitative influence on the superconductivity 
properties of the system (12ii] . 

Let us now calculate the entanglement in the BCS 
state. We will focus on the entanglement in momen- 
tum space, i.e., entanglement between all electrons hav- 
ing quantum numbers a and their time-reversed electrons 
described by —a. We use the logarithmic ne gati vity, a 
proper measure of entanglement, to this aim UJ. The 
negativity of a quantum state p, M{p), is the sum of the 
absolute value of the negative eigenvalues of the partial 
transpose . The logarithmic negativity is then given 
bylAA(p) = log(l + 2AA(p))p. 

From equation (6) is easily seen that the negativity 
between the electrons described by a and —a is given 
by Afa,-a = UaVa- Using the additivity property of the 
logarithmic negativity |13l | we then find that the entan- 
glement between all the electrons described by a and the 
time-reversed electrons described by —a can be written 
as 



EM{BCS) = Y,^og[l 



[(e„-/i)2 + A2]i/2 



(7) 



One finds that entanglement is a monotonic increasing 
function of the order parameters, having, therefore, quan- 
titative influence in the superconductivity properties of 
the system. Actually, the logarithmic negativity by itself 
could be used as an order parameter. 

The rj pairing model: Despite the great success of the 
BCS theory in describing superconductivity at very low 
temperatures, it fails in explaining the superconductiv- 
ity of some materials at high temperatures. Actually, the 
mechanism which allows the existence of superconductiv- 
ity at temperatures as high as 160K is not yet completely 
understood. We now analyze a particular mechanism for 
Tc superconductivity, the r] pairing of electrons Q . The 
biggest difference between this and the BCS model is that 
in the former Coopers pairs are formed from electrons at 
the same site, whereas in the later electrons forming a 
Cooper pair have an average flnite separation distance. 



Some of the entanglement properties of the r] pairing 
model were studied in Refs. 0, 0, 0, 0|. It was 
particularly shown that again in this model the existence 
of entanglement is a necessary condition for supercon- 
ductivity JJs J/Q]. We will expand this result, establish- 
ing relations between the amount of entanglement and 
the order parameter of the system in both the flnite and 
thermodynamical regimes. 

Following Ref. [13, consider a set of sites, where each 
one can be occupied by fermions having spin down and 
up. Let cj g be the fermionic creation operator, where the 
indices i and s identify the «-th site and the spin orienta- 
tion, respectively. The operators c satisfy the following 
commutation relations 



(8) 



The operator that creates a coherent superpositions of 
Cooper pairs in each of the sites, 77^, is given by 



i,] id' 



(9) 



where n is the number of sites of the system. The op- 
erators can be applied several times, where in each 
one a new superposition is created. Nonetheless, due to 
the Pauli exclusion principle, the number of applications 
cannot exceeds the number of sites. The state of a sys- 
tem in which k coherent pairs were created is given by 



\k, n — k) 



(10) 



Note that in this representation each site is described ef- 
fectively by one qubit, whose value stands for an empty 
site and 1 for a occupied one . 

The most important characteristic of the 77 states is 
the existence of off diagonal long range order (ODLRE), 
which implies the main superconductor properties, such 
as the Meissner effect and the flux quantization 17]. The 
ODLRE is defined by a non-zero value of the off diagonal 
elements of the reduced density matrices of two sites, 
when the distance between these two becomes arbitrarily 
long: 



lim 



a. 



(11) 



where a is a constant independent of n in the thermo- 
dynamical limit. The number a is exactly the order pa- 
rameter of the high-Tc superconductivity . 

The reduced density matrix aij for the sites i and j of 
the state |n, n — fc) is 

ay = a|00)(00| + fe|ll)(ll| + 2c|^+)(V'+|, (12) 

where = (|01) + |10))/^/2, 

fc(fc- 1) 



n{n — 1) ' 



(n — /s)(n — fc — 1) 

h^- T 13 

n{n — 1) 
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and 



2k{n- k) 
n{n — 1) 



(14) 



Since the state is symmetric, this density matrix is the 
same no matter how far the two sites are from each other. 
The order parameter is thus easily found to be a = c. In 
the next paragraphs we associate it with the entangle- 
ment of the system. In this model we study the entan- 
glement between the occupation number of each site. 

Consider first the case of finite n. That is the case, for 
instance, of superconductivity in nano-structures |23| . In 
Ref. 16] it was shown that, in this case, the two sites 
density matrices are entangled for every n and fc, and 
that the order parameter c is just the maximum fidelity 
of teleportation under local operations and classical com- 
munication i22|llBl- Notably, we see that in the finite 
case the number c quantifies both the superconductivity 
features of the system and its usefulness as a quantum 
channel! 

For general superconducting materials we must con- 
sider the thermodynamical regime, where n, fc — *■ oo. In 
this scenario the analysis becomes more complex, as the 
states tJij become separable no matter what are the ratio 



between n and k |l6l Il8|. From equation (11) we find 
that the order parameter is non-zero iff 



lim — 

n,fc — ^oo il 



(15) 



for some real number r. As noted by Vedral [T^ [l7j |. 
although there is no bipartite entanglement between two 
sites in this regime, the system still has multipartite en- 
tanglement [231 . In order to quantify it, we will use 
the logarithm geometric measure of entanglement LRq, 
which is given by the logarithm of the overlap between 
the state and its nearest separable state in the norm 2 
sense j23|. This measure was calculated for the states 
\n,n— k) and is given by [U 



LRa{n, k) = log 



k\{n~k)\ 



(?)' 



-k 



Using k — rn and the Stirling approximation ln{nl) 
nlnn — n, we have 



LEG{n,k) = log 



(l + r) 




One might note that when — s- cx), the amount of entan- 
glement diverges. However, if we use instead the density 
of entanglement, i.e., the average entanglement per site, 

cLe ~ LEc/n^ we find 



dE — lim 



LEain.k) 



-r log r — (1 — r) log(l — r) 



From equation (15) one finds in addition that a = 
2r(l — r). Therefore, it holds the following direct re- 
lation between the amount of multipartite entanglement 
and the order parameter 



dE 



log 



1 + 



log 



1 + VT 



(16) 



Interestingly, the relation between dE and a is same re- 
lation of the entanglement of formation and the concur- 
rence j22| , which implies that also in this model the order 
parameter quantifies entanglement. 

The Bose-Hubbard Model: The simplest non-trivial 
model for interacting bosons in a periodic potential is 
described by the Bose-Hubbard (BH) Hamiltonian [23j. 
It contains most of the physics of strongly interacting 
bosons, i.e., the competition between kinetic and inter- 
action energy. The BH Hamiltonian is given by [23l | 



Hbh 



-J 



<Ti,m> n 



(K-M)a,\an, 
(17) 



where is the annihilation operator for an atom at site 
TO. The first term in equation (17), proportional to J, is 
the tunneling matrix element between nearest-neighbors. 
The parameter U is proportional the repulsion intensity 
of two atoms at the same site. 

A zero temperature, the physics of the Bose-Hubbard 
model can be divided in to extreme regimes. The first is 
the one where J is much lower than U and the system is 
described by a Mott insulator. The other is the regime 
dominated by the kinetic energy, where J is much larger 
than U and the system presents superfluidity properties 

M 

The ground state of the Bose-Hubbard Hamiltonian 
cannot be found analytically for every U and J. There- 
fore, we will restrict our analysis to the limiting cases 
f7/J — > oo and U/J 0. Consider a system with N 
atoms. When U/J 0, its ground state is a deeply 
superfluid state given by |2^ 



IV'sf) 




(18) 



where M is the number of wells of the system. In the 
other limit, where the system is described by a Mott 
insulator, the ground state is just a separable state given 

by 



IM 



M 



(19) 



where \gi) are variational local states. The order param- 
eter usually considered in this model is the expectation 
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value of the annihilation operator am In the ther- 

modynamical limit N, M ^ oo, with 



Thus we can replace the sum in equation (24) by the 
integral 



N,M^oo V M 



(20) 



it can be shown that in the MI phase (a„i) = 0, whereas 
in the superfluidity phase (am) = r 

Consider now the bipartition where one party is formed 
by the Fock space associated with the operator and 
the other by the Fock space associated with all other op- 
erators an, n ^ m. The state |?/'sf) can then be written 
as 



\tpSF) 




(21) 



with & = Y J2n^m '^n- Equation (21) can expressed 
in terms of its Schmidt decomposition 

1^^^^) = E {(^P^^^ -p)'^\k)ajN~k),, (22) 

where p = 1/Af. Let us calculate the entanglement in 
this partition, which quantifies the quantum correlations 
existent between the number of atoms in site m and the 
number of atoms on all the other sites. Using the nega- 
tivity |l3] as a measure of entanglement, 



AA(^sf) = 



N-k 



(23) 



From the central limit theorem we have that in the limit 
^ oo. 



g2]vp(i-p)_ (24) 



(25) 

Therefore, in this limit, we have again a direct relation 
between entanglement and the value of the order param- 
eter: 



lim Afi-ipsp) = lim 2— 1 , , 



N 



1 



1/2 



1 



V2r--. 



In conclusion, we have shown that there exist a direct 
relation between entanglement and the order parame- 
ter in three important physical systems. This connec- 
tion has important consequences both in the theory of 
phase transitions and in quantum information science. 
On one hand, quantum order parameters gain an inter- 
esting physical interpretation, which clarifies the intrinsic 
quantum character of QFT. On the other hand, the rich 
literature about the properties of QOPs, including evolu- 
tion equations and statistical properties, can be used in 
the study of entanglement and the viability of quantum 
information processing in those systems. 

Note added: After this work was completed, I learned 
that a connection between single-site entanglement and 
the energy-gap in the BCS model was obtained by Shi in 
Ref. 0- In this same paper an interesting relation be- 
tween entanglement and the proper fractional part of the 
filling factor is also derived, which extends the discussion 
of the present paper also to the quantum Hall state. 
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